REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 

Davis  Highway.  Suite  1 204.  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington.  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

January  2008  USARIEM  Technical  Report 

4.  TITLE  AND  SUBTITLE 

Corrections  to  the  Shapiro  Equation  Used  to  Predict  Sweating  and  Water 

Requirements 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

R.R.  Gonzalez,  S.N.  Cheuvront,  D.A.  Goodman,  L.A.  Blanchard,  L.G.  Berglund, 
M.N.  Sawka 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Thermal  and  Mountain  Medicine  Division 

U.S.  Army  Research  Institute  of  Environmental  Medicine 

Natick,  MA  01760-5007 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Same  as  #7  above 

10. SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  words) 

This  project  tested  the  original  equation  (OSE)  in  which  sweating  rate,  msw  (g  *  mA-2  *  hA-l)  =  27.9  *  Ereq  * 

EmaxA-0.455,  Shapiro  et  al.  (EJAP  48:  83,1982).  OSE  predicts  msw  over  wide  heat  stress  conditions  (HS),  clothing  systems 
including  chemical  protective  (CPS),  and  metabolic  activities.  Ereq  is  required  sweat  loss  calculated  from  heat  balance  and 
(Emax)  is  the  maximum  evaporative  power  of  a  thermal  environment.  OSE  is  suitable  over  a  limited  HS  up  to  2-h  and  work 
rates  <  450W.  Another  objective  was  to  develop  a  new  algorithm  or  correct  the  OSE  so  that  reliable  fluid  replacement 
guidelines  can  be  constructed  for  work  times  >  2-h.  USARIEM  studies  revealed  that  OSE  in  a  Heat  Str  ess  Decision  Aide 
(HSDA)  over-predicts  actual  msw  over  wide  HS,  work  rates,  and  times.  Overpredictions  in  msw  often  lead  to  over-hydration 
problems.  Data  from  101  volunteers  (80  men  and  21  women)  were  analyzed  including  experiments  at  various  activity  levels 
over  wide  HS  (Ta  =  15°C  -  46°C,  Pw=  2  -  33  Torr  and  (V)  =  0.4  -  2.5  m*sA-l)  and  various  CPS  with  body  armor. 

Studies  included  4  separate  chamber  studies,  1  USARIEM  field  study  &  1  study  outside  the  U.S.  Data  were  analyzed  using 
fuzzy  piecewise  regression.  OSE  predicted  high  values  (p<  0.003)  compared  to  observed  msw.  We  recommend  HSDA 
incorporate  the  algorithm:  msw  =  147+  1.527*(Ereq)  -  0.87*(Emax),  which  accounts  for  effects  of  heavy  work,  CPS,  body 
armor,  8h  work,  and  is  not  gender  specific.  A  correction  to  OSE  is:  msw  =  147*exp  (0.0012*OSE).  These  equations  require 
testing  over  wider  effective  radiant  loads  in  the  field  (Solar  Load). 

14.  SUBJECT  TERMS 

Sweat  prediction,  physiological  modeling,  water  requirements 

15.  NUMBER  OF  PAGES 

44 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

Unclassified 

NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18  298-102 


USAPPC  VI  .00 


DISCLAIMER 


Approved  for  public  release:  distribution  is  unlimited. 

The  opinions  or  assertions  contained  herein  are  the  private  views  of  the  author(s) 
and  are  not  to  be  construed  as  official  or  reflecting  the  views  of  the  Army  or  the 
Department  of  Defense. 

The  investigators  have  adhered  to  the  policies  for  protection  of  human  subjects 
as  prescribed  in  Army  Regulation  70-25,  and  the  research  was  conducted  in 
adherence  with  the  provisions  of  32  CFR  Part  219. 

Human  subjects  participated  in  these  studies  after  giving  their  free  and  informed 
voluntary  consent.  Investigators  adhered  to  AR  70-25  and  USAMRMC 
Regulation  70-25  on  the  use  of  volunteers  in  research. 

Any  citations  of  commercial  organizations  and  trade  names  in  this  report  do  not 
constitute  an  official  Department  of  the  Army  endorsement  of  approval  of  the 
products  or  services  of  these  organizations. 


USARIEM  TECHNICAL  REPORT  T07-## 


CORRECTIONS  TO  THE  SHAPIRO  EQUATION  USED  TO  PREDICT 
SWEATING  AND  WATER  REQUIREMENTS 


Richard  R.  Gonzalez1 
Samuel  N.  Cheuvront2 
Daniel  A.  Goodman2 
Laurie  A.  Blanchard3 
Larry  G.  Berglund3 
Michael  N.  Sawka2 


biology  Department,  New  Mexico  State  University,  Las  Cruces,  NM 
2Thermal  and  Mountain  Medicine  Division 
3Biophysics  and  Biomedical  Modeling  Division 


January  2008 


U.S.  Army  Research  Institute  of  Environmental  Medicine 
Natick,  MA  01760-5007 


TABLE  OF  CONTENTS 


Section  Page 

List  of  Figures . iv 

List  of  Tables . v 

Background . vi 

Acknowledgments . vii 

Executive  Summary . 1 

Introduction  . 2 

Methods  . 4 

Results  . 11 

Discussion  . 24 

Conclusions . 29 

Recommendations . 29 

References . 30 

Appendix  A . 34 

Appendix  B . 36 


LIST  OF  FIGURES 


Figure  Page 


Oxygen  uptake  (VO2)  from  the  data  set 

1A  12 

Sweating  rate  from  all  experiments  and  individuals 
IB  14 

Residual  line  plot  using  original  Shapiro  equation  (OSE)  vs. 

2  observed  data  for  sweat  loss  (g*m'2*h'1)  15 

Output  from  OSE,  exponential  correction  to  OSE  and  observed 

3  sweating  rate  plotted  against  piecewise  regression  output.  17 

Negative  exponential  smoothing  denoted  by  solid  lines  to  the 

data 

Sweat  loss  (g*m'2*h'1)  output  (means  ±95%  confidence 

4  intervals)  from  the  various  equations  and  observed  data  plotted  18 

vs.  the  given  trials 

Predicted  sweat  rate  responses  plotted  vs.  observed  sweat  rate 

5  for  all  trial  data  19 

Output  from  the  various  prediction  equations  and  observed 
6A  sweat  rate  plotted  for  each  subject  in  a  Canada  NBC  study  20 

Output  from  each  prediction  equation  plotted  vs.  observed 
6B  sweating  rate  21 

Output  from  each  of  the  prediction  equations  and  observed 
7A  sweat  loss  plotted  for  each  individual  walking  in  MOPP  clothing  22 

Log  plot  of  the  outputs  from  the  prediction  equations  vs. 

7B  observed  sweat  loss  23 

Original  data  from  Shapiro  estimated  by  a  linear  regression  of 
8  Ereq/dmsw  vs.  Emax  (r2=0.929)  25 


IV 


LIST  OF  TABLES 


Table  Page 

1  Key  descriptive  and  physiological  data  for  2-h  experiments  in  6 

Protocol  H03-14  used  in  developing  the  present  algorithm 

2  Key  descriptive  and  physiological  data  for  8-h  experiments  in  7 

Protocol  H03-14  used  in  developing  the  present  algorithm 

3  Heat  production  (M,  W*m'2)  and  observed  sweat  loss  (OSL,g*m'  13 
2*h_1).  Values  are  least  squares  means  between  men  (M)  and 
women  (F)  for  combined  data 


V 


BACKGROUND 


A  person’s  fluid  balance  and  thermoregulatory  equilibrium  are  essential  for  peak 
performance.  Dehydration  to  levels  exceeding  2%  normal  body  mass  and/or  elevation 
of  body  temperature  above  38.5°C  are  associated  with  reduced  physiologic  function  and 
cognitive  performance,  and  these  impairments  increase  with  increasing  levels  of  body 
water  loss,  particularly  in  hot  environments.  Verification  of  valid  and  accurate  fluid 
replacement  algorithms  for  individuals  (including,  but  not  limited  to  firefighters,  police 
personnel,  athletes,  and  Warfighters)  is  a  critical  step  in  the  development  of  “smart,” 
personal  guidelines  and  computer  devices  that  can  be  implemented  to  optimize 
performance  and  prevent  heat  injuries.  The  Institute  of  Medicine  (IOM  2005)  has 
identified  the  “development  of  capabilities  to  predict  hourly  and  daily  water  requirements 
based  on  metabolic  rate,  climatic  conditions,  and  clothing”  as  a  research  priority. 

Adequate  hydration  and  core  temperature  equilibrium  are  essential  for  the 
prevention  of  heat-related  illnesses  and  for  sustaining  peak  physical  and  cognitive 
performance.  Hydration  is  particularly  important  during  training  and  field  operations, 
especially  under  hot  climatic  conditions.  Evaluation  of  the  current  operational  water 
requirement  guidelines  using  a  systematic  physiologic  approach  will  enhance  a  better 
decision-making  process  regarding  activity  in  hostile  environment  conditions  faced  by 
both  military  and  civilian  populations. 

During  regular  military  activity  there  are  about  120  heat-stroke  injuries  each  year, 
which  are  associated  with  a  $1  OM/y  cost  (3,31 );  see  also  http://asma.army.mil).  In  the 
military  population,  dehydration  is  a  co-morbidity  factor  in  20%  or  more  reported  heat 
casualties  (3).  Validations  of  fluid  requirements  and  physiological  strain  assessment 
improve  safety  and  performance,  facilitate  mission  planning,  optimize  hydration  and 
logistical  efficiency,  and  help  personnel  overall  in  conduction  operational  missions. 
Furthermore,  evaluation  of  the  current  guidelines,  obtained  from  studies  performed  on  a 
wider  database,  should  address  broader  physiologic  issues  pertinent  to  both  military 
and  civilian  operational  needs  where  there  is  currently  a  gap  in  operational  knowledge. 

In  addition,  this  study  effort  should  provide  a  better  decision-making  process  regarding 
activity  in  hostile  environment  conditions,  which  eventually  will  help  in  decreasing  the 
expenses  associated  with  physical  and  heat  stroke  injuries.  Algorithms  that  are  easily 
migrated  into  current  rational  and  empirical  models  addressing  correct  fluid 
requirements  and  physiological  strain  assessment  can  improve  the  safety  and 
performance  of  Warfighters  and  optimize  hydration  and  logistical  efficiency. 
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EXECUTIVE  SUMMARY 


The  objective  of  this  project  was  to  cross  validate  the  original  equation, 
OSE,  in  which  sweating  rate  (g*m'2*h'1)  =  27.9  •  Ereq  •  (Emax)’0  45  ,  which  was 
developed  by  Shapiro  et  al.  (33)  to  predict  rate  of  sweat  loss  over  wide 
environmental  conditions,  clothing  systems,  and  metabolic  activities.  Ereq  is 
required  sweat  loss  calculated  from  heat  balance,  and  Emax  is  the  maximum 
evaporative  power  of  a  thermal  environment.  Within  the  limits  of  the  data,  the 
equation  has  been  shown  to  be  a  valid  estimator  of  sweating  rate  for  a  variety  of 
heat  stress  exposures  up  to  2  h  and  work  rates  limited  to  less  than  450  W.  The 
second  objective  was  to  develop  a  new  prediction  algorithm  or  correct  the 
previous  one  so  that  reliable  fluid  replacement  guidelines  using  such  equations 
can  be  constructed  in  the  future  for  more  extended  work  times  greater  than  just  2 
h,  the  exposure  time  present  in  the  original  study  database.  The  need  to  develop 
a  new  prediction  equation  stems  from  results  of  recent  studies  at  USARIEM  that 
revealed  that  the  original  equation  embedded  in  a  Heat  Strain  Decision  Aide 
(HSDA)  computer  model  over  predicts  actual  sweating  responses  over  wide 
environmental  extremes,  work  rates,  and  work  periods.  Overpredictions  of 
sweating  rate  (and  the  required  fluid  intake  to  fully  replace  the  expected  sweat 
loss  during  extremes  of  heat  loss)  can  lead  to  over-hydration  problems.  A 
USARIEM  database  was  secured  consisting  of  101  volunteer  subjects  (80  men 
and  21  women)  who  completed  experiments  at  various  activity  levels  over  wide 
environmental  ambient  conditions.  Ta  ranged  from  15°C  to  46°C,  ambient  water 
vapor  pressures  varied  from  Pw=  2  to  33  Torr,  and  air  movements  (V)  were  from 
0.4  to  2.5  nrs'1.  Subjects  wore  various  military  clothing  systems  including 
chemical  protective  clothing  and  body  armor.  Raw  data  were  obtained  from  4 
separate  chamber  studies  and  1  field  study  conducted  at  USARIEM,  and  1 
laboratory  study  conducted  by  the  Defence  R&D  Canada  (DRDC),  Toronto.  Each 
element  of  the  comprehensive  heat  balance  equation  was  analyzed.  Data  were 
analyzed  using  fuzzy  piecewise  linear  and  nonlinear  regression  analyses  to 
establish  appropriate  change  points  in  sweat  loss  per  time  points.  It  was 
established  that  the  original  Shapiro  algorithm  tested  in  this  study  (101  subjects, 
longer  work  durations  from  4-8  h,  and  a  variety  of  clothing  systems)  predicts 
markedly  high  values  in  sweating  rates.  The  most  important  finding  of  the  current 
study  and  recommendation  is  to  substitute  or  modify  the  current  HSDA  program 
with  a  corrected  algorithm:  sweat  loss  (g*m'2*h'1)  =  147+1 .527*(Ereq)  - 
0.87»(Emax).  This  equation  takes  into  account  effects  of  heavy  work  and  clothing 
factors,  body  armor,  longer  exposure  times  (8  h),  and  is  not  gender  specific. 
Alternatively,  a  correction  to  the  original  equation  can  be  used  as  a  simple 
replacement:  sweat  loss  =  147»exp  (0.001 2*OSE).  These  equations  require 
testing  over  wider  effective  radiant  loads  in  the  field  (effect  of  Solar  Load)  using  a 
larger  database. 
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INTRODUCTION 


Maintaining  heat  balance  by  regulating  core  temperature  within  tolerable 
limits  by  having  proper  hydration  levels  over  complex  thermal  environments  is 
essential  for  peak  human  performance.  Dehydration  levels  exceeding  2%  of 
normal  body  mass  coupled  by  elevations  in  core  temperature  above  38°C  are 
associated  with  reduced  physical  and  cognitive  performance  (6,7,16,30,31,34), 
particularly  in  hot  environments.  Likewise,  core  temperatures  below  36°C  can 
also  lead  to  hypothermia  and  cold  injury. 

The  U.S.  military  currently  uses  two  prediction  models,  based  on  two 
different  concepts,  to  predict  sweating  rates  and  core  temperatures  of  Soldiers 
during  operational  stress.  The  U.S.  Army  Research  Institute  of  Environmental 
Medicine  (USARIEM)  Heat  Strain  Model  (HSDA)  is  an  empirical  model  that 
includes  equations  to  predict  sweating  rate,  work/rest  cycles,  and  maximal 
working  times  during  different  levels  of  metabolic  rate  and  exercise  intensity 
(13,15,27).  These  algorithms  have  been  used  to  prepare  military  guidance 
regulations  for  water  needs  and  work/rest  cycles  during  training  and  deployment 
(e.g.,  FM  10-52;  FM  21-10;  TBMED-507).  Recently,  these  algorithms  were 
embedded  in  a  digital  mini  heat  strain  monitor  (HSM),  which  is  a  rugged,  pocket- 
sized,  advisory  heat  strain  device  that  also  outputs  integrated  WBGT  and, 
therefore,  may  be  applied  in  the  field  or  for  various  training  maneuvers. 

Predictive  equations  for  implementing  work/rest  disciplines  with  various  clothing 
systems,  environments,  or  workload  sequences  in  the  model  are  based  on  a 
series  of  functions  that  cascade  into  final  output  values  of  work/rest  cycles,  water 
requirements,  and  maximum  endurance  times  for  a  given  environmental  activity 
and  clothing  system  scenario. 

The  second  thermoregulatory-cardiovascular  model,  called  SCENARIO 
(1 1,18),  is  a  rational  Soldier-physiologic  model  designed  to  simulate  the  time 
course  of  heat  strain  observed  during  military,  industrial,  and  athletic  settings.  In 
this  model,  the  human  body  is  modeled  as  a  single  cylinder  containing  six 
compartments:  [1]  a  central  core  representing  the  heart,  lungs,  and  splanchnic 
regions,  [2]  a  muscle  layer,  [3]  a  subcutaneous  fat  layer,  [4]  a  vascular  skin  layer, 
[5]  a  superficial  vascular  skin  layer,  and  [6]  a  central  blood  compartment.  Rates 
of  metabolic  heat  production  in  the  core,  fat,  and  vascular  skin  compartments  are 
assumed  to  be  fixed  percentages  of  total  resting  metabolism.  Heat  production  of 
the  muscle  layer  is  variable,  depending  on  total  energy  expenditure  and  external 
work  performed.  It  is  assumed  that  the  blood  and  skin  compartments  produce  an 
insignificant  amount  of  heat.  The  SCENARIO  model  (11)  can  be  easily  migrated 
into  a  robust  military  operational  computer  simulation  software  product  for  use  in 
planning  and  mission  evaluation  of  Soldier's  individual  variability. 

In  1972,  Givoni  and  Goldman  (15)  developed  an  empirical  approach  to 
predict  core  temperature  response.  They  inferred  that  for  any  given  combination 
of  metabolic  rate,  environment,  and  clothing,  a  theoretically  determined 
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equilibrium  maximal  point  of  core  and  skin  temperatures  would  be  generated, 
and  unified  equations  could  be  constructed  to  predict  that  set  limit.  A  series  of 
predictive  equations  were  developed  and  were  proved  useful  in  describing 
human  heat  exchange  while  subjects  wore  a  variety  of  clothing  systems  that 
matched  with  core  temperature  pattern  response  during  rest,  work,  and  recovery 
in  the  heat. 

The  predictive  equations  for  implementing  work/rest  cycles  with  various 
clothing  systems,  climates,  or  workload  sequences  in  the  original  Givoni- 
Goldman  model  were  based  on  a  series  of  theoretical  regression  equations  that 
simulated  the  thermoregulatory  function  (27).  The  main  equation  in  the  model 
was  one  that  established  the  difference  in  core  temperature  expected  at 
theoretical  equilibrium  fora  given  environment,  work  rate,  and  clothing  system. 
The  core  temperature  at  every  minute  was  dependent  on  the  initial  core 
temperature,  the  equilibrium  core  temperature,  and  the  time  in  the  period  from 
the  end  of  a  proceeding  delay  time  during  work  or  rest.  A  significant  number  of 
training  guidelines  for  work/rest  cycles  in  the  U.S.  Army,  the  Israeli  Defense 
Forces,  the  Technical  Cooperation  Program  (TTCP)  countries,  and  many  other 
military  training  guidelines  were  based  on  these  modified  equations  (13). 

Crucial  to  the  above  model’s  operational  application,  and  embedded  in  its 
model  computer  code,  is  the  Shapiro  sweating  rate  algorithm  currently  used  in 
the  HSDA,  a  modern  sequel  of  the  original  Givoni-Goldman  model.  The  latter 
model  was  developed  from  limited  laboratory  experiments  on  men  only  for 
energy  expenditures  ranging  from  approximately  75  W  (rest)  up  to  475  W 
(moderate  intensity  work  for  a  dismounted  Soldier)  over  a  range  of  environmental 
conditions  (20°-54°C  and  10%-94%  relative  humidity  [RH])  while  wearing  shorts 
and  T-shirt,  outdated  military  fatigues,  and  CB  garments.  It  is  also  only  applicable 
to  predict  sweating  responses  for  men  (33).  The  derived  equation  is  shown 
below: 


Sweating  rate  (g*rrf2‘h'1)  =  27.9  •  Ereq  •  (Emax)'0455  (Eq.  1 ) 

where  Ereq  is  the  evaporation  required  to  maintain  heat  balance  at  any  given  core 
temperature  and  Emaxis  the  maximal  evaporative  capacity  of  the  environment. 
This  equation  has  been  used  to  predict  water  requirements,  assuming  the  fluid 
intake  (L/h)  replaces  the  expected  water  lost  by  sweating  in  a  fully  heat- 
acclimated  person  (sweating  rate  x  body  surface  area  (BSA)  x  1 0'2,  L/h). 

A  caveat  to  the  use  of  model  analyses  that  overpredict  sweat  losses  is 
that  excessive  and  persistent  overdrinking  can  lead  to  a  relatively  rare  but 
potentially  lethal  condition  known  as  symptomatic  hyponatremia  or  water 
intoxication  (6,10,23,26).  Current  guidelines  found  in  various  military  operational 
manuals  have  been  shown  to  overpredict  requirements.  Use  of  inaccurate  sweat 
loss  predictions  can  lead  to  incorrect  hydration  needs  (forced  or  permissible) 
during  extreme  heat  and  exercise  challenges  (4,22,23).  An  outbreak  of  exercise- 


3 


induced  hyponatremia  cases  in  the  U.S.  military  (10)  led  to  a  revision  of  the 
recommended  hourly  fluid  consumption  guidelines  during  work  in  hot  weather 
(22). 


While  the  array  of  exercise  and  environmental  conditions  would  appear  to 
have  covered  adequately  a  broad  range  of  climatic  conditions  using  the  2-h 
Shapiro  equation  to  predict  water  requirements,  only  resting  exposures  were 
included  at  air  temperatures  below  35°C.  Additionally,  the  clothing  evaluations 
were  done  using  static  manikin  (still  air)  evaluations  and  do  not  cover  modern 
clothing  materials  now  being  tested  (employing  dynamic  manikins)  and  used  by 
military  developers  for  the  modern  Warfighter.  Because  accurate  prediction  of 
water  requirements  is  essential  for  more  than  just  mission  requirements  in  hot 
climates,  it  is  imperative  that  physiologic  algorithms  be  accurately  validated  to 
predict  sweating  rates  over  wider  environments,  including  both  temperate  and 
cool  climatic  zones  over  which  Warfighters  conduct  operations. 

The  purposes  of  this  study  were  the  following:  [1]  to  assess  the  accuracy 
of  the  original  Shapiro  sweat  loss  prediction  equation  (established  solely  using  a 
2-h  time  frame)  using  a  wider  database  that  includes  extended  time  periods 
during  prolonged  exercise  (8  h),  heavy  work  rates,  and  contemporary  clothing 
systems  including  body  armor;  [2]  to  seek  to  develop  new  algorithms  that  can 
predict  sweating  requirements  that  are  more  accurate  over  extended  periods, 
and  [3]  to  provide,  if  possible,  a  correction  to  the  original  Shapiro  equation  (OSE) 
that  is  applicable  and  easily  migrated  into  the  current  HSDA  for  use  during  more 
extended  periods  of  exercise  in  cool  and  warm  environments.  The  hypotheses, 
based  on  current  research  and  literature  results  (4,22,23),  were  that  OSE 
overestimates  sweating  rates  over  extended  periods  of  exercise  and,  therefore, 
leads  to  excessive  estimates  of  water  requirements  for  particular  operational 
needs.  A  new  algorithm  based  on  a  more  extensive  database  than  the  original 
would  allow  a  more  finite  estimation  of  water  requirements  over  wider 
environments  and  exercise  intensities.  Additionally,  a  correction  to  the  old 
equation  should  allow  a  reasonable  prediction  of  the  water  requirements  without 
scrapping  a  25-year  history  of  operational  doctrine  related  to  water  requirements, 
and  should  prevent  unnecessary  water  loading.  Such  an  improved  algorithm 
would  be  a  useful  feature  for  implementation  in  present  HSDA  water  requirement 
predictions  decisive  in  future  military  and  civilian  operational  applications. 


METHODS 


DATA  ANALYSIS 

The  database  consisted  of  101  volunteer  subjects  (80  men  and  21 
women)  who  completed  experiments  at  various  activity  levels  over  wide 
environmental  ambient  conditions.  Ta  ranged  from  15°C  to  46°C,  ambient  water 
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vapor  pressures  varied  from  Pw=  2  to  33  Torr,  and  air  movements  (V)  were  from 
0.4  to  2.5  nrs'1.  Subjects  wore  various  military  clothing  systems.  Each  protocol 
was  approved  by  the  appropriate  institutional  review  boards,  and  all  volunteers 
were  informed  both  verbally  and  in  writing  of  the  objectives  and  procedures  of  the 
respective  study.  No  identifications  of  a  given  volunteer’s  personal  records  were 
present  in  the  current  spreadsheet  database. 

Raw  data  were  obtained  from  4  separate  chamber  studies  and  1  Field 
study  conducted  at  USARIEM,  and  1  laboratory  study  conducted  by  the  Defence 
R&D  Canada  (DRDC),  Toronto  (13).  These  studies  are  described  in  the  following 
sections: 

(a)  USARIEM  protocol  H98-04  (Montain,  PI).  Data  were  examined  from  19 
previously  heat  acclimatized  Soldiers  (13  men  and  6  women)  who  completed  all 
experiments  dressed  in  hot  weather  Battle  Dress  Uniforms  (BDUs,  clo  value 
1.08;  im/clo=  0.49).  Extensive  study  details  can  be  found  in  Montain,  et  al.(22).  In 
brief,  individuals  completed  12  randomize  exercise-heat  stress  trials  in  which 
they  walked  at  3  exercise  intensities  either  at  250,  425  and  600  W  in  three  humid 
environments  (Ta=28°C/Pw=21 .3  Torr;  Ta=32°C/Pw=26.8  Torr  and 
Ta=36°C/Pw=33.4  Torr).  In  other  studies  they  walked  at  425  W  in  three  dry 
environments  (  Ta=36°C/Pw=1 1  Torr;  Ta=41°C/Pw=  14.6  Torr;  and 
Ta=46°C/Pw=19  Torr).  Dry  heat  stress  trials  were  completed  following  a  humid 
test  condition.  Appropriate  work/rest  cycles  for  each  exercise  task  were  initially 
determined  using  the  HSDA  (27)  over  a  2-h  total  exposure.  Rectal  temperatures, 
heart  rate,  and  oxygen  uptake  from  indirect  calorimetry  (2-min  Douglas  bag 
collection)  were  measured.  Missing  mean  skin  temperature  (Tsk)  data,  required  to 
calculate  pertinent  heat  balance  equation  parameters,  was  estimated  using 
Saltin’s  (28)  equation:  Tsk=0.215*Ta+26.6  (±0.5  SEE).  Whole  body  sweating 
rates  were  calculated  from  changes  in  pre-  and  post-exercise  body  weights  after 
correction  for  clothing  weights,  water  intake,  urine  production,  and  estimates  of 
metabolic  water  loss  (mres)  and  respiratory  heat  loss  (Eres)  (see  Appendices  A  & 
B). 


(b)  USARIEM  protocol  H03-14  (Montain,  PI).  Details  can  be  found  in 
Cheuvront  et  al.  (4).  Thirty-nine  healthy  individuals  participated  in  this  study.  The 
clothing  ensemble  was  the  U.S.  Army  woodland  BDU  with  field  cap,  sleeves 
down  (clo=1 .08,  im/clo=0.49  at  wind  speed  1  m/s),  and  athletic  shoes.  Test 
sessions  lasted  either  2-  or  8-h.  Twenty-one  volunteers  (16  men  and  5  women) 
participated  in  the  2-h  experiments,  and  physical  characteristics  for  this  group  are 
shown  in  Table  1 .  Eighteen  different  volunteers  (17  men  and  1  woman) 
completed  the  8-h  experiments,  and  their  characteristics  are  shown  in  Table  2. 
Tables  1  and  2  also  describe  the  7  different  levels  of  environmental  stress,  work: 
rest  cycles,  VO2  at  each  work  intensity,  and  sweat  loss  observed.  In  the  2-h 
study  trials,  volunteers  were  not  heat  acclimated,  while  in  8-h  experiments, 
volunteers  were  previously  heat  acclimated. 
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Table  1.  Key  descriptive  and  physiological  data  for  2-h  experiments  in  Protocol  H03-14  used  in  developing  the  present 
algorithm.  _ _ _ _ _ _ _ _ 


Trial 

Ta 

(°C) 

Pw 

(Torr) 

Work:Rest 

Cycles 

#(min) 

N 

Males 

Females 

BSA 

(m2) 

Body 

Weight 

(kg) 

Tsk 

(°C) 

V02 

(L/min) 

SR  (L/h) 

A1 

15 

6.4 

2x  (50:10) 

12M 

1.94±0.15 

77.6±8.8 

29.3  ±  1.0 

0.94  ±0.14 

0.135  ±0.079 

4F 

1.70±0.03 

64.3±4.2 

28.2±0.82 

0.79±0.17 

0.188±0.086 

A 

15 

6.4 

2x  (50:10) 

13M 

1.96±0.14 

79.7±9.3 

28.7  ±  0.8 

1.49  ±0.17 

0.305  ±0  .136 

2F 

1.68 

60.4 

28.7 

1.3 

0.319 

B 

15 

6.4 

2x  (50:10) 

13M 

1.96±0.14 

79.7±9.3 

29.1±1.0 

1.99±0.21 

0.472±0.170 

2F 

1.68 

6.4 

27.5 

1.57 

0.424 

C 

20 

8.8 

2x  (50:10) 

14M 

1.96±0.14 

80.7 

31.1  ±0.8 

178  ±  15 

0.220  ±  .089 

IF 

1.63 

56.9 

31.34 

0.81 

0.165 

D 

20 

8.8 

2x  (50:10) 

15M 

1.98±0.13 

80.9±8.9 

30.2±1.08 

1.45±0.19 

0.410±0.177 

4F 

1.67±0.04 

62.6±5.7 

30.4±0.47 

1.14±0.13 

0.346±0.084 

E 

20 

8.8 

2x  (50:10) 

15M 

1.97±0.14 

80.5±9.4 

30.3±1.15 

2.03±0.26 

0.625±0.199 

2F 

1.69 

63.5 

29.6 

1.49 

.498 

F 

25 

11.9 

2x  (50:10) 

11M 

1.96±0.1 1 

80.2±8.9 

31.8±0.55 

1.03±0.16 

0.321±0.102 

2F 

1.68 

60.35 

31.1 

0.84 

0.421 

G 

25 

11.9 

2x  (50:10) 

9M 

1.95±.12 

80.0±9.5 

31.5±0.92 

1.47±0.20 

0.479±0.171 

IF 

1.73 

63.8 

30.9 

1.26 

0.515 

H 

25 

11.9 

2x  (50:10) 

10M 

1.95±0.1 1 

80.5±8.9 

32.1±0.62 

1.99±0.27 

0.755±0.237 

IF 

1.63 

56.9 

31.6 

1.54 

0.558 

1 

30 

15.9 

2x  (50:10) 

10M 

1.95±0.1 1 

80.5±8.9 

33.5±0.75 

1.98±0.26 

0.935±0.317 

IF 

1.63 

56.9 

33.1 

1.45 

0.702 

Legend:  n  (number  of  subjects  and  gender  completing  experiment),  Ta  (air  temperature),  Pw  (ambient  water  vapor  pressure),  Tsk  (mean  skin 
temperature),  SR  (observed  sweating  rate).  All  data  are  means  ±  SD  except  when  n  <3.  Wind  speed  for  all  experiments  was  1  m/s. 
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Table  2.  Key  descriptive  and  physiological  data  for  8-h  experiments  in  Protocol  H03-14  used  in  developing  the  present 
algorithm. 


Trial 

Ta 

(°C) 

Pw 

(Torr) 

Work:  Res 
t  Cycles 
#(min) 

N 

Males 

Females 

BSA 

(m2) 

Body 

Weight 

(kg) 

Tsk 

(°C) 

V02 

(L/min) 

SR  (L/h) 

J 

40 

22.1 

6x  (60:20) 

12M 

2.03±0.14 

84.7±12.4 

35.9  ±  0.61 

1.11  ±0.21 

0.667  ±0.114 

IF 

1.63 

55.3 

36.2 

0.775 

0.404 

K 

35 

12.7 

6x  (60:10) 

15M 

1.97±0.12 

80.2±9.8 

34.4  ±  0.5 

1.39  ±0.21 

0.569  ±0  .062 

IF 

1.62 

54.8 

34.12 

0.92 

0.405 

L 

35 

12.7 

6x  (60:10) 

15M 

1.99±0.15 

81.7±12.9 

33.9±0.5 

1.05±0.15 

0.452±0.058 

M 

27 

10.7 

6x  (60:10) 

14M 

1.98±0.15 

81.7±12.9 

32.3  ±  0.6 

1.44  ±0.23 

0.406  ±  0.096 

IF 

1.63 

55.7 

31.83 

0.96 

0.263 

N 

27 

10.7 

6x  (60:10) 

12M 

1.98±0.16 

81.1±13.9 

32.7±0.64 

1.05±0.14 

0.269±0.05 

IF 

1.63 

55.6 

33.0 

0.770 

0.153 

0 

20 

8.8 

6x  (60:10) 

12M 

1.97±0.16 

80.9±14 

30.6±1.15 

1.39±0.22 

0.229±0.087 

IF 

1.63 

55.1 

30.3 

0.847 

0.162 

Legend:  n  (number  of  subjects  and  gender  completing  experiment),  Ta  (air  temperature),  Pw  (ambient  water  vapor  pressure),  Tsk  (mean  skin 
temperature),  SR  (observed  sweating  rate).  All  data  are  means  ±  SD  except  when  n  <3.  All  subjects  previously  heat  acclimated.  Wind  speed  for 
all  experiments  was  1  m/s. _ 
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(c)  USARIEM  protocol  H05-12  (Cheuvront,  PI).  Details  can  be  found  in 
Cheuvront  (5).  Thirteen  men  completed  3  experimental  trials  in  a  hot,  dry 
environment  (Ta=  35°C,  Pw=12.7  Torr,  Im/s  wind  speed).  The  U.S.  Army  BDU 
was  worn  in  all  3  trials,  either  alone  (trial  BDU),  combined  with  Interceptor  Body 
Armor  (trial  IBA),  or  combined  with  Interceptor  Body  Armor  and  spacer  vest  (trial 
SP).  Four  hours  of  intermittent  treadmill  walking  was  performed  (50  min  walking: 
10  min  rest  intervals).  In  the  BDU  trial,  walking  speed  was  1 .56  m/s  with  a  grade 
of  3%.  In  the  IBA  and  SP  trials,  the  grade  was  reduced  to  2%  to  compensate  for 
the  added  armor  weight,  which  allowed  examination  of  the  clothing  contributions 
to  heat  strain  independent  of  added  load  carriage. 

(d)  H06-14  (Cadarette,  PI).  One  female  and  five  male  volunteers  did 
continuous  exercise  (treadmill  speed  at  80.2  m/min;Vo2  ranges  0.86  to  1 .24 
L/min)  for  2  h  while  dressed  in  BDU  +  IBA,  as  in  (c)  (IBA:  clo  =  1 .35,  im/clo  =  0.27 
at  wind  speed  1  m/s).  Environmental  conditions  were  Ta=30°/Pw=16  Torr; 
Ta=35°C/Pw=32  Torr;  and  Ta=40°C/  Pw=1 1  Torr.  All  procedures  were  as 
outlined  in  (c)  above. 

Clothing 

All  clothing  characteristics  (insulation,  clo;  vapor  resistance,  im/clo)  were 
measured  on  a  sweating  thermal  manikin  for  each  clothing  configuration  tested 
(BDU,  IBA,  SP).  In  the  BDU  trial,  the  BDU  was  worn  with  field  cap,  sleeves  down, 
and  athletic  shoes,  rather  than  field  boots,  to  reduce  blister  formation  (BDU:  clo  = 
1.12,  im/clo  =  0.44  at  wind  speed  of  1  m/s).  In  the  IBA  trial,  subjects  also  wore  the 
Interceptor  Body  Armor  vest  to  include  front  and  rear  ballistic  protective  inserts 
(throat  and  groin  protection  excluded).  The  outer  vest  is  made  of  a  fine  Kevlar 
weave,  and  the  protective  plates  of  boron  carbide  ceramic  with  spectra  shield 
backing.  The  total  weight  of  the  vest  as  used  was  7.5  kg,  and  it  covered  -25%  of 
the  total  BSA  (IBA:  clo  =  1 .35,  im/clo  =  0.27  at  wind  speed  1  m/sec).  In  the  SP 
trial,  subjects  also  wore  a  1-cm  thick  vest  of  proprietary  knit  fabric  in  between  the 
Interceptor  Body  Armor  and  uniform.  The  spacer  vest  is  designed  to  produce  an 
air  channel  that  can  theoretically  increase  the  potential  for  ventilation  and 
evaporative  cooling  of  the  torso  (SP:  clo  =  1 .28,  im/clo  =  0.32  at  wind  speed  1 
m/sec).  The  weight  of  the  spacer  vest  was  nominal  (0.2  kg). 

Procedures 


At  the  start  and  conclusion  of  each  trial,  nude  body  mass  was  measured 
on  an  electronic  precision  balance  scale  (Toledo  1 D 1 ,  Worthington,  OH; 
accuracy  ±  20g).  To  minimize  differences  in  hydration  state  between  trials,  250 
ml  water  was  given  1  h  before  starting  exercise.  Body  mass  measurements  that 
fluctuated  by  less  than  1%  of  the  10-day  mean  were  considered  normal.  A 
similar  hydration  state  from  trial  to  trial  was  assumed.  Additional  fluid  was  given 
on  the  morning  of  a  test  only  if  body  mass  deviated  downward  by  more  than  1% 
of  the  10-day  mean  (5). 
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Heart  rate  (HR)  (Polar  a3,  Polar  Electro,  Inc.,  Woodbury,  NY)  and  core 
(intestinal)  body  temperature  (Tc)  (Jonah™  core  body  temperature  capsule,  Mini 
Mitter  Company,  Inc,  Bend,  OR)  were  measured  continuously  and  recorded  at 
10-min  intervals. 

Independent  Cross  Validation  Data  Sets: 

Ft.  Bliss  Field  Study.  A  separate  archival  raw  data  set  was  obtained 
from  a  previous  field  study  conducted  at  Ft.  Bliss,  TX,  by  USARIEM  and  used  to 
test  the  algorithms.  Details  can  by  found  in  Santee  et  al.  (29).  In  brief,  8  males 
(average  body  mass=80.5±1 5.2  kg  SD;  BSA  =  1 ,97±0. 1 8  m2  SD)  walked  at  a 
pace  of  2  mph  for  12  miles  (24  min  continuous  exercise  with  a  6-min  break)  on  a 
calibrated  track  for  6  h.  Subjects  carried  a  22-kg  pack,  and  average  heat 
production  was  maintained  at  194.5±20.4  W*m'2.  Clothing  systems  were  MOPP  0 
(BDU:  clo=  1 .34/  im/clo=0.31  at  Im/s),  MOPP  1  (clo=1 .97/  im/clo=0.17  at  1  m/s), 
and  MOPP  4  (clo  =  2.44/  im/clo=0.12).  Volunteers  did  night  and  day  walks  in  the 
various  clothing  configurations.  Core  temperature  was  measured  with  a 
thermistor  inserted  10  cm  into  the  rectum,  and  3-site  skin  temperatures  were 
measured  with  heat  flow  discs;  all  data  were  collected  using  a  Grant  squirrel 
temperature  recording  device.  Daytime  ambient  conditions  were  Ta=23°C / 

Pw=1 .6  Torr,  with  a  mixed  solar  load  estimated  by  the  Effective  Radiant  Field  at 
500W*m'2,  resulting  in  an  operative  temperature  (T0)  =  49.5°C  (9). 

DRDC,  Toronto  Data  Set.  Raw  data  sets  from  a  completed  laboratory 
study,  as  part  of  a  report  to  The  Technical  Cooperation  Program  (TTCP)  (13), 
were  used  to  test  the  various  prediction  equations.  Data  from  13  males  and  9 
females  (follicular  phases  of  their  menstrual  cycle)  were  used  to  compare  the 
various  equations.  Physical  characteristics  (±SD)  were  males,  body  mass, 
82.7±12.5  kg;  BSA,  2.01±0.16;  females,  body  mass,  60.4±8.9;  BSA,1 ,66±0.15. 
The  ambient  conditions  were  Ta  =  40°C/  Pw=  1 1  Torr  at  V=0.4m*s'1.  Subjects 
exercised  for  2  h  (4  work:rest  cycles  of  15:15  min)  or  before  their  rectal 
temperature  reached  a  peak  value  of  39.5°C,  or  heart  rate  became  elevated  no 
higher  than  1 80  beatsTnin'1  for  3  min.  Average  M  (±SD)  of  the  male  groups  was 
203. 4±  24.1  W*m'2  and  for  the  female  group,  average  M  was  187.5±17.4  W*m'2. 
Rectal  temperature  was  measured  from  a  thermistor  inserted  15  cm  beyond  the 
anal  sphincter  and  mean  skin  temperature  from  a  12  point  area  weighted 
average  recorded  on  strategic  skin  area  sites.  Subjects  were  dressed  in  Canada 
NBC  clothing  system  previously  evaluated  by  USARIEM  manikin  procedures  and 
found  to  have  a  clo  value  of  1 .88/  im/clo=  0.18  at  V=1m*s'1  (14). 
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HEAT  TRANSFER  ANALYSIS  PROCEDURES 


Each  element  of  the  comprehensive  heat  balance  equation  (Appendices  A 
&  B)  was  determined  from  the  raw  data  and  concatenated  to  link  together  in  a 
unified  spreadsheet  (Microsoft  Excel®)  for  later  model  algorithm  analysis  using 
various  statistical  modules  (STATISTICA®,  Version  7,  Tulsa,  OK)  and  nascent 
proprietary  mathematical  code. 

STATISTICAL  PROCEDURES 

Quasi-Newton  Method 


All  concatenated  data  were  analyzed  using  fuzzy  piecewise  linear  and 
nonlinear  regression  analyses  (35,36)  to  establish  appropriate  change  points  in 
sweat  loss  per  time  points  in  the  data  set,  coded  by  trial,  sex,  and  individual 
subject  number,  and  derive  intercepts  for  independent  parameters  derived  from 
the  heat  balance  equation  (Ereq,  Emax  )  for  the  data  set.  A  Quasi-Newton  method 
was  employed  to  derive  regression  parameters.  In  this  method,  the  slope  of  a 
function  at  a  particular  point  is  first  computed  as  the  first-order  derivative  of  the 
function  (at  that  point).  The  "slope  of  the  slope"  becomes  the  second-order 
derivative,  which  reveals  how  fast  the  slope  is  changing  at  the  respective  point, 
and  in  which  direction.  The  quasi-Newton  method  will,  at  each  step,  evaluate  the 
function  at  different  points  in  order  to  estimate  the  first-order  derivatives  and 
second-order  derivatives.  The  analysis  technique  then  uses  this  information  to 
follow  a  path  towards  the  minimum  of  the  loss  function  (35).  The  fuzzy  piecewise 
is  more  robust  than  conventional  methods  and  not  sensitive  to  outliers  or 
irregular  data,  as  was  found  in  the  present  data  set  (see  Figures  1  and  1A).  The 
technique,  as  constructed  for  this  study,  is  also  suitable  for  long-term  time  series 
predictions  of  specific  variables  (35).  Comparison  of  the  original  Shapiro  equation 
predicting  sweat  loss  and  water  requirements  was  initially  done  against  observed 
data  to  obtain  residual  analyses  to  ascertain  how  much  the  Shapiro  equation 
deviated  from  the  observed  data.  Next,  the  new  fuzzy  piecewise  equation  (35) 
was  compared  against  the  original  Shapiro  equation  and  the  observed  raw  data 
secured  for  each  separate  trial.  Corrections  to  the  original  Shapiro  equation  were 
derived  by  independent  piecewise  regression  analyses  incorporating  an  iterative 
approach  to  obtain  the  most  optimum  equation  (exponential,  log  fit,  etc.)  and  test 
the  significance  of  the  derived  regression  coefficients  (Wald  statistic)  that  fit  the 
database.  The  conditions  were  that  the  new  equation  prediction  of  sweat  loss  did 
not  deviate  from  the  observed  data  or  the  independently  determined  fuzzy 
piecewise  equation  by  more  than  ±0.28  L/h  (roughly  ±150  g*m'2*h'1). 

Data  are  expressed  as  means  ±SD  or  ±SE,  or  as  means  ±95%  confidence 
interval  (Cl).  The  differences  in  observed  sweat  loss,  heat  production,  and  the 
output  from  the  various  prediction  equations  were  analyzed  by  a  Factorial 
ANOVA  design  to  include  main-effects  and  interactions  for  categorical  predictors 
(gender,  all  Trials).  Both  univariate  (using  a  given  single  continuous  dependent 
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variable)  and  multivariate  (multiple  continuous  dependent  variables)  designs 
were  analyzed.  If  a  significant  F  value  was  found  for  a  given  dependent  variable, 
the  Bonferroni  adjustment  procedure  was  used  as  a  post  hoc  approach  to  locate 
critical  differences  at  P<0.01  and  considered  statistically  significant;  correlations 
with  a  probability  value  greater  than  .01  (including  those  with  p-values  between 
.01  and  .05)  were  considered  non-significant. 

Finally,  an  independent  cross  validation  analysis  of  the  fuzzy  piecewise 
equation  was  executed  against  two  independent  archival  data  sets:  a  field  study 
in  which  a  group  of  men  walked  in  various  levels  of  MOPP  clothing  (0,  I,  and  IV) 
(29)  and  using  data  in  a  lab  study  conducted  at  DRDC,  Toronto  (13)  previously 
discussed  in  the  Methods  section. 

RESULTS 

Figure  1 A  shows  the  VO2  obtained  from  each  of  the  subjects  during  the 
various  trials.  Along  with  the  other  individuals’  parameters  (Tsk,  clo,  im/clo)  and 
environmental  parameters,  heat  production  from  the  raw  VO2  was  determined  by 
the  comprehensive  heat  balance  to  obtain  Ereq  (Appendices  A  and  B). 

Figure  1 B  shows  observed  sweating  rates  for  each  trial  in  the  data  set 
depending  on  exercise  intensity,  ambient  condition,  and  clothing  system  worn. 
The  observed  sweating  rates  were  compared  with  the  original  Shapiro  equation 
(Eq.  1 )  to  determine  residuals. 
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Figure  1A.  Oxygen  uptake  (V02)  from  the  data  set. 
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The  Ereq  and  the  Emax  obtained  by  solution  of  the  heat  balance  equation  for  each 
individual  and  specific  environmental  parameter  were  used  to  derive  a  new 
equation  using  fuzzy  piecewise  regression,  as  described  in  the  Methods  section. 

Table  3  shows  the  calculated  mean  heat  production  (M,  W*m'2)  and 
observed  sweat  loss  (OSL,  g*m'2*h'1)  ±SE  for  each  trial  (Methods)  separated  by 
gender  for  the  USARIEM  data  set.  Following  ANOVA  tests,  which  might  reveal 
whether  there  was  any  significant  effect  of  gender,  a  Bonferroni  post  hoc  (20) 
analysis  was  done  on  the  data  set.  There  was  significant  difference  in  OSL  only 
between  men  and  women  within  the  moderate  intensity  exercise  trials  of  Protocol 
H98-04,  combining  all  dry  and  wet  environments. 
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Table  3.  Heat  production  (M,  W*m'2)  and  observed  sweat  loss  (OSUgnrT^h'1).  Values  are 
least  squares  means  between  males  (M)  and  females  (F)  for  combined  data. 


Cell 

Trial 

Sex 

M 

M 

(— SE) 

OSL 

OSL 

(-SE) 

N 

1 

H98-04-L 

M 

163.04 

8.87 

295.21 

18.35 

32 

2 

H98-04-L 

F 

148.60 

12.18 

208.31 

25.17 

17 

3 

H98-04-M 

M 

224.00 

5.84 

407.13 

12.06** 

74 

4 

H98-04-M 

F 

252.24 

8.37 

332.59 

17.29** 

36 

5 

H98-04-H 

M 

329.18 

8.25 

443.96 

17.06 

37 

6 

H98-04-H 

F 

357.61 

11.83 

356.88 

24.46 

18 

7 

H03-14-B 

M 

252.42 

3.53 

224.72 

7.30 

202 

8 

H03-14-B 

F 

231.82 

10.04 

210.19 

20.75 

25 

9 

H05-12-B 

M 

284.66 

13.92 

353.21 

28.78 

13 

10a 

H05-12-B 

F 

NA 

NA 

NA 

NA 

11 

H05-12-I 

M 

287.79 

13.92 

407.64 

28.78 

13 

12a 

H05-12-I 

F 

13 

H05-12-S 

M 

289.52 

13.92 

420.47 

28.78 

13 

14a 

H05-12-S 

F 

15 

H06-14-I 

M 

213.07 

12.96 

362.12 

26.79 

15 

16 

H06-14-I 

F 

193.80 

28.98 

483.49 

59.91 

3 

**  P<0.01 ,  Bonferroni  post-hoc  test  within  trial  cells;  all  others  NS.  SE=  standard  error;  N=  number  of  values  within 
each  group.  aNot  estimable  due  to  lack  of  or  small  sample  size. 
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Figure  1 B  shows  the  observed  sweating  rates  during  all  experiments  that 
were  used  to  compare  against  the  original  Shapiro  equation  and  analyzed  to 
develop  a  new  prediction  algorithm. 


Figure  IB.  Sweating  rate  from  all  experiments  and  individuals. 

USARIEM  Dataset 
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Individual  Experimental  Data  Values 


Noticeable  are  the  extensive  high  points  and  depressions  reflecting  the  variable 
sweating  outputs  of  the  individuals  during  the  mixed  work  intensities  and 
environmental  conditions. 
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Figure  2  shows  the  results  of  the  residual  analysis  comparing  the  output 
from  the  original  Shapiro  equation  against  the  observed  sweating  shown  in 
Figure  1A. 


Figure  2.  Residual  line  plot  using  output  from  the  original  Shapiro  equation  vs. 
observed  data  for  sweat  loss  (g^m'2^'1). 

Residual  Line  Plot  (Old  Shapiro  Eq  vs  Obs  Sw  Loss) 


It  is  clear  from  Figure  2’s  plot  of  the  residual  values  (comparison  of  a  given 
observed  sweat  loss  value  minus  the  predicted  value  for  each  data  point)  that  the 
original  Shapiro  equation  overpredicted  sweat  loss  by  as  much  as  100%  from  the 
various  trials,  especially  during  heavy  work  intensities  (H98-04),  and  under¬ 
predicted  by  some  80%  during  cooler  trials  (Ta=15°C).  However,  there  is  a  range 
in  the  data  in  which  the  model  predicts  within  ±20%  of  observed  data  particularly 
during  easy  work  activities  and  mild  heat  stress  conditions. 

The  data  set  was  next  examined  to  develop  optimum  regression 
parameters  that  would  satisfy  all  experimental  trials  sufficiently  (within  an  SEE  of 
±0.28  L/min  [~±150  g^m'2^'1]  and  coefficient  of  determination,  R>0.8).  Since  the 
data  were  collected  at  different  times  and  contained  disparate  conditions 
(environmental  stress,  gender,  acclimation  states,  work  intensities  from  heavy  to 
easy  and  various  clothing  systems),  a  non-linear  fuzzy  piecewise  regression 
employing  a  Quasi-Newton  solution  (see  Methods)  gave  the  best  resolution 
(35,36).  ^ 

The  reduced  fuzzy  piecewise  regression  algorithm  (fPW)  developed  for 
the  combined  data  set  from  a  solution  of  Ereq  and  Emax  is: 
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Sweating  rate  (g*m'2*h'1)  =  147  +  1 .527*Ereq  -  0.87*Emax  (Eq.  2) 

resulting  in  an  R  =  0.89,  which  explained  78.94%  of  the  variance  in  the  disparate 
data  set.  For  this  analysis,  Ereq  and  Emax  were  determined  from  the  individual 
heat  balance  equation  and  transformed  by  division  of  W*m'2/(0.68  W'h'^g'1). 

Since  it  was  determined  that  the  output  from  the  original  Shapiro  equation 
(OSE,  Eq.  1)  would  probably  not  unify  the  data  set  and  predict  adequately 
(Figure  2),  a  similar  iterative  approach  used  to  obtain  Eq.  2  was  run  to  correct  the 
OSE  using  individual  analysis  of  Ereq  and  Emax  in  the  data  set.  An  exponential 
correction  to  Eq.  1  was  successfully  obtained  after  attempting  various  statistical 
algorithm  solutions.  This  solution  produced  the  following  correction  equation  to 
predict  sweating  rates: 

Correction  to  sweating  rate  (g*m'2*h'1)  =  147»exp^0,0012*OSE^  (Eq.  3) 

where  OSE  in  Eq.  3  is  the  uncorrected  output  from  the  original  Shapiro  equation 
(33): 


A  msw  =  27.9*Ereq*Emax'0'455  (Eq.  1) 

Figure  3  shows  a  comparison  of  the  OSE,  the  exponential  correction  to 
OSE,  and  the  observed  sweating  rates  from  the  present  data  set  plotted  against 
output  values  calculated  with  the  new  piecewise  regression  sweat  prediction 
equation.  The  solid  lines  indicate  negative  (decreasing)  exponential  smoothing 
curves  secured  by  a  polynomial  regression  algorithm.  In  the  latter  regression,  the 
weights  that  determine  the  influence  of  individual  data  points  of  the  sweat  loss  on 
consecutive  segments  of  a  curve  are  calculated  using  a  negative  exponential 
function.  Using  this  function,  unnoticeable  patterns  of  data  response  can  be 
revealed.  For  instance,  the  OSE  response  curve  is  unremarkable  within  the 
range  of  20-300  g*m'2*IT1  and  is  probably  the  optimum  range  that  this  equation 
should  be  applied;  however,  it  begins  to  deviate  from  the  corrected  OSE 
prediction  equation  and  the  observed  sweat  loss  data  curves  at  about  320  g*m" 
^•h'1.  There  is  a  noticeable  “break”  in  the  OSE  curve  occurring  around  500  g*m' 
2*h_1,  which  continues  to  a  high  ceiling  point  that  is  almost  double  the  value 
denoted  by  the  piecewise  regression  prediction  output.  Both  the  exponential 
correction  to  OSE  output  and  the  observed  data  match  the  piecewise  regression 
output  precisely  up  to  a  value  of  500  g*m'2*h'1.  Past  this  point,  the  exponential 
correction  to  the  OSE  equation  begins  to  deviate  upwards  while  the  observed 
data  curve  exhibits  a  downward  trend.  The  latter  phenomenon  may  be  perhaps 
owing  to  low  sweating  responses  during  low  work  intensity  and/or  cool 
environmental  trial  runs  (Ta=15°C).  Linear  regression  analyses  (shown  in  the  box 
insert  to  Figure  3)  estimated  by  plotting  the  outputs  of  the  3  dependent  variables 
(OSE,  corrected  OSE,  and  observed  sweat  loss  data)  against  the  output  from 
piecewise  equation  show  high  rvalues  (0.96,  0.88,  and  0.78,  respectively). 
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However,  these  r-values  cannot  quantitatively  verify  the  individual  segmentation 
in  response. 

Figure  3.  Output  from  OSE,  exponential  correction  to  OSE  and  observed 
sweating  rate  plotted  against  piecewise  regression  output.  Negative  exponential 
smoothing  is  denoted  by  solid  lines  in  the  data. 


USARIEM  Data  Exponential  Smoothing 
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Figure  4.  Sweat  loss  (g*m'2*h'1)  output  (means  ±95%  confidence  intervals)  from 
the  various  equations  and  observed  data  plotted  vs.  the  given  trials.  Values  are 
staggered  for  clarity.  Double  asterisk  (**)  for  the  OSE  prediction  shows  significant 
differences  between  trials  and  right  bracket  (])  between  other  curves  at  p<0.003 
(Bonferroni).  There  is  a  non-significant  difference  between  the  other  prediction 
equations  and  observed  sweat  loss  for  each  trial. 


USARIEMData 
Mean±0.95  Conf.  Interval 


H98-04-M  H03-14-B  H05-12-I  H06-14-I 

Trial 


Figure  4  above  shows  a  comparison  of  the  output  from  the  various 
prediction  equations  along  with  the  observed  sweat  rates  lumped  for  each  trial 
run.  The  OSE  showed  the  widest  variability  compared  to  the  observed  data  and 
was  markedly  elevated  for  the  moderate  and  heavy  work  intensity  experiments. 
However,  the  corrected  prediction  equation  tracked  output  response  remarkably 
well  in  comparison  to  the  observed  data  and  piecewise  prediction  equation. 
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Figure  5.  Predicted  sweat  rate  responses  plotted  vs.  observed  sweat  rate  for  all 
trial  data. 


USARIEM  Data 
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Exp  Pred  Correction  OSE  =  136. 8476+0. 6027*x;  0.95  Conf.lnt. 


Observed  Sweat  Rate  (g-irf^h'1) 


Figure  5  shows  the  output  from  the  various  prediction  equations  plotted 
against  observed  data.  Noticeable  is  the  fact  that  the  output  from  the  original 
Shapiro  equation  overpredicts  by  almost  66%  above  500  g*m'2»h'1,  and  below 
that  range  the  scatter  is  widespread.  Using  the  corrected  equation  compares  well 
with  the  piecewise  regression  equation,  but  the  scatter  is  still  prevalent. 
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Cross  Validation  Analysis  from  Two  Separate  Studies 


Figure  6A.  Output  from  the  various  prediction  equations  and  observed  sweat  rate 
plotted  for  each  subject  in  a  Canada  NBC  study. 
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To  see  how  well  the  various  prediction  equations  track  individual  response 
data,  a  comparison  test  was  done  using  separate  study  results:  one  using  field 
data  and  the  other  using  a  laboratory  study  conducted  outside  USARIEM.  The 
outputs  from  the  various  prediction  equations,  as  well  as  the  observed  sweat 
loss,  are  shown  in  Figure  6A. 

The  OSE  predicted  sweat  loss  response  higher  than  observed  data 
consistently  across  subjects  (p<0.003),  but  so  did  output  generated  from  the 
piecewise  regression  equation  (by  about  50%  less).  It  is  not  clear,  at  this  time, 
why  this  occurred.  One  reason  may  be  that  the  experimental  runs  during  this  trial 
were  conducted  at  chamber  wind  speed  of  0.4  m/s,  while  the  dynamic  im/clo 
evaluated  for  the  clothing  system  had  been  at  1  m/s.  The  evaporative  heat 
transfer  coefficient  strongly  coupled  with  skin  wettedness  (8)  (cohe,  Wnn'^Torr'1) 
and  other  heat  balance  parameters  in  the  original  piecewise  regression 
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coefficients  that  are  used  to  calculate  Ereq  and  Emax  are  co-mingled  with  free  and 
forced  convection  and  possibly  not  sensitive  to  low  wind  speed.  Therefore,  latent 
heat  transfer  within  layers  of  chemical  protective  clothing  is  calculated  at  a  higher 
value.  Interestingly,  the  corrected  equation  to  the  OSE  matched  the  observed 
sweating  rate  responses  between  subjects. 

Figure  6B.  Output  from  each  prediction  equation  plotted  vs.  observed  sweating 
rate. 


Canada  NBC  Data 


Observed  Sweating  Rate  (g*nT2*h'1) 


Figure  6B  shows  this  fact  more  clearly  when  the  output  from  each  of  the  3 
prediction  equations  is  regressed  against  the  observed  sweating  rates.  The 
corrected  equation  matches  the  observed  data.  However,  in  the  linear  plot 
analysis,  the  regression  coefficient  (1.427)  of  the  piecewise  regression  equation 
is  significantly  higher  (p<0.003)  than  the  one  evident  from  the  corrected  OSE 
(0.912).  This  indicates  that  the  latent  heat  loss  may  be  elevated  falsely 
particularly  between  the  interstices  of  military  chemical  protective  clothing  layers. 
Additional  research  is  required  to  investigate  this  property. 
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Figure  7A.  Output  from  each  of  the  prediction  equations  and  observed  sweat  loss 
plotted  for  each  individual  walking  in  MOPP  clothing. 
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The  outputs  from  the  various  prediction  equations  and  observed  sweat 
loss  during  a  field  study  test  (29)  are  shown  in  Figure  7A.  In  general,  the 
piecewise  regression  and  corrected  OSE  equations  predicted  individual  sweating 
rates  adequately  (within  the  ±150  gnrf2*h'1  criterion  deviation)  except  for  a  few 
individuals.  The  OSE,  however,  consistently  predicted  too-high  values. 

This  response  is  more  clearly  evident  in  Figure  7B  in  which  the  above  data 
are  transformed  into  logarithmic  form  and  plotted  in  linear  coordinates.  The  OSE 
regression  coefficient  (1 17.6)  is  significantly  higher  than  the  other  two  regression 
coefficients  (p<0.001 ).  The  regression  coefficients  from  the  outputs  using  the 
piecewise  regression  equation  and  the  corrected  OSE  are  not  significantly 
different  and  each  follows  the  observed  data  more  closely  than  the  OSE. 
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Figure  7B.  Log  plot  of  the  outputs  from  the  prediction  equations  vs.  observed 
sweat  loss. 


Piecewise  Eq  =  216.6012+51 ,8649*log10(x) 
Original  Shapiro  Eq  =  213.9182+1 17.5599*log10(x) 
Shapiro  Correction  =  172.891 +40. 2289*log10(x) 
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DISCUSSION 


The  objective  of  this  project  was  to  cross  validate  (using  new  data)  the 
original  equation  that  was  developed  by  Shapiro  et  al.  (33)  to  predict  rate  of 
sweat  loss  over  wide  environmental  conditions,  clothing  systems,  and  metabolic 
activities,  which  did  not  include  women  responses  in  its  formulation.  The 
rationale  to  develop  a  predictive  equation  to  gauge  sweat  loss  (and,  thereby, 
water  requirements)  was  a  unique  concept  at  the  time.  This  equation  has  been 
used  to  predict  water  requirements  over  wide  thermal  environments  by 
knowledge  of  only  two  key  variables,  Ereq  and  Emax,  which  directly  or  indirectly 
integrate  the  effects  of  the  internal  factors  (M,  skin  and  core  temperature)  and 
external  factors  (clothing,  operative  temperature,  wind,  and  humidity).  Within  the 
limits  of  the  data,  the  equation  has  been  shown  to  be  a  valid  estimator  of 
sweating  rate  for  a  variety  of  heat  stress  exposures  of  up  to  2  h  and  work  rates 
limited  to  less  than  450  W. 

The  next  objective  was  to  develop  a  new  prediction  algorithm,  or  correct 
the  previous  one  so  that  reliable  fluid  replacement  guidelines  using  such 
equations  can  be  constructed  in  the  future  for  more  extended  work  times  greater 
than  just  2  h,  the  exposure  time  present  in  the  original  study  database.  The  need 
to  develop  a  new  prediction  equation  stems  from  results  of  recent  studies 
(4,5,22;  Xu,  personal  communication,  2007)  that  revealed  that  the  original 
Shapiro  equation  embedded  in  a  Heat  Strain  Decision  Aide  (HSDA)  computer 
model  tended  to  overpredict  actual  sweating  responses  over  wide  environmental 
extremes,  work  rates,  and  work  periods.  The  consequences  are  that 
overpredictions  of  sweating  rate  (and  the  required  predicted  fluid  intake  to  fully 
replace  the  expected  sweat  loss  especially  during  extremes  of  heat  loss)  can 
lead  to  over-hydration  problems  (10,23,26). 

To  accomplish  these  objectives,  a  wider  database  was  necessary  that 
includes  heavy  exercise,  modern  clothing  systems  including  body  armor,  and 
incorporated  sweating  responses  from  both  men  and  women.  The  present  study 
includes  such  data  (over  100  subjects). 

A  critical  shortcoming  was  found  in  the  data  that  produced  the  sweat  loss 
prediction  equation  (Eq.  1)  selected  from  the  original  reference  data  (33).  This 
scrutiny  revealed  that  the  equation  was  developed  by  a  relationship  of  Ereq/  Amsw 
plotted  against  Emax  using  averaged  data  from  three  series  of  experiments  in 
which  Amsw  is  sweat  loss  (g*m'2*h'1).  This  regression  equation 
(Y=0.0537(Emax)0  455)  was  then  solved  as  Amsw  =  Ereq/  0.0537(Emax)0  455  to  obtain 
18.6  (the  reciprocal  of  0.0537)  x  Ereq/  (Emax)0  455  and  converted  by  division  of  the 
latent  heat  of  vaporization  of  sweat  (A,,  J/g)  to  obtain  the  final  equation  present  in 
HSDA  model.  A  linear  approximation  would  have  satisfied  the  averaged  data, 


24 


and  it  is  not  certain  why  this  was  not  pursued  in  the  data  set,  as  evident  in  Figure 
8. 


Figure  8.  Original  data  from  Shapiro  estimated  by  a  linear  regression  of  Ereq 
/drrisw  vs.  Emax  (r2=0.929).  In  addition,  the  solution  of  the  above  equation  for 
dmsw  is  714Ereq/(Emax  +190),  W/m2,  or  1 050Ereq/(Emax  +279)  g*m'2*h'1. 


Shapiro  Data,  (ref  33) 

Erq/dmsw  =  0.2656+0.0014*x;  0.95  Pred.lnt. 


Emax  (W/m2) 


Another  key  shortcoming  of  such  analysis  evident  from  the  original  report 
(see  Figure  4,  Shapiro  et  al.,  (33))  is  the  use  of  average  data  solely  to  obtain  a 
final  equation,  which  inherently  biases  the  residual  sum  of  squares  (lessens  its 
variance)  in  a  regression  analysis  shown  in  Figure  8.  Also,  the  7,  unit  used  was 
incorrect,  which  Wenger  (37)  showed  is  2,426  J/g  (0.68  W*h/g),  thereby  resulting 
in  the  sweat  loss  prediction  equation  being: 

Sweating  rate  (g-rn'^h'1),  Amsw  =  27.4  ■  Ereq  ■  (Emax)'0455  (Eq.  4) 

Additionally,  the  assumption  is  incorrect  that  Ereq/  Amsw  is  not  a  correlate  of  M 
(metabolic  heat  production),  when  in  fact  Ereq  is  based  on  the  solution  of  the 
comprehensive  heat  balance  equation  comprising  M-W-(R+C)-  S.  Another 
limitation  in  the  original  sweat  loss  equation,  that  is  now  embedded  in  the  HSDA, 
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is  the  capitation  of  mean  skin  temperature  to  a  value  of  36.5°C,  rather  than 
allowing  skin  temperature  to  “float”  at  the  specific  level  for  a  given  work  intensity 
and  environmental  challenge  where  steady-state  heat  balance  can  be  obtained. 
This  skin  temperature  (and  skin  saturation  vapor  pressure,  Ps,Sk)  increases  dry 
heat  (R+C)  and  latent  heat  losses  (elevating  the  Ps,skto  ambient  vapor  pressure 
gradient  through  contiguous  clothing  layers)  and  artificially  raises  the  predicted 
sweating  rates  at  a  given  time  point,  possibly  initiating  an  overprediction  in  total 
sweat  loss. 

In  the  present  study,  output  from  the  original  Shapiro  prediction  equation 
consistently  overestimated  measured  rate  of  sweat  loss  during  high  intensity 
exercise  up  to  a  maximum  in  every  trial  from  both  the  present  database  and  the 
archival  data.  Maximal  sweating  rate  is  generally  limited  to  about  667  g*m'2*h'1,  a 
value  close  to  1 .27  L/h  for  a  person  with  a  BSA  of  1 .9  m2  (38,39).  Above  that 
value,  a  person  generally  cannot  reach  values  necessary  to  achieve  steady-state 
heat  balance. 

A  reevaluation  of  the  original  Shapiro  algorithm  using  the  present 
extensive  database  (101  subjects,  longer  work  durations  from  4-8  h,  and  a 
variety  of  clothing  systems)  established  that  the  original  equation  for  Amsw  does 
indeed  statistically  predict  markedly  high  values  in  sweating  rates.  Therefore, 
water  requirements  using  the  equation  would  also  predict  too-high  fluid 
requirements  and  concur  with  recent  studies  (4,22)  showing  similar  over¬ 
estimations  using  that  equation. 

The  most  important  accomplishment  of  the  current  analyses  was  the 
development  of  an  improved  equation: 

msw  (g'm'2^'1)  =  147  +1 ,527*(Ereq)  -0.87»(Emax)  (Eq.  2) 

based  on  fuzzy  piecewise  regression  analysis  that  incorporates  the  combined 
effects  of  wider  metabolic  intensities,  body  armor  and  clothing  systems,  and  is 
essentially  applicable  for  both  men  and  women  working  for  longer  time  periods 
than  2  h  if  proper  work:rest  cycles  are  interposed.  The  equation  still  incorporates 
attributes  of  the  original  Shapiro  equation:  one  related  to  solution  of  the  heat 
balance  equation  (Ereq)  and  the  other  associated  with  clothing  worn  and 
environmental  impact  (Emax). 

Additionally,  a  correction  to  the  original  equation  modeled  from  similar 
concepts  was  constructed  in  which: 

msw  (g*m'2*IT1)  =  147*exp(0.0012*OSE)  (Eq.  5) 

where  OSE  is  output  from  the  original  equation  that  may  be  employed  to  predict 
water  requirements  to  replace  fluid  lost  from  sweating  within  ±0.125  mL/h.  The 
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limits  of  the  above  predictive  equations  comprise  the  original  Ereq  and  Emax  limits 
(50<Ereq<360  W*m'2and  20<Emax  <525  W*m'2,  respectively),  but  extends  these 
limits  so  the  algorithm  is  applicable  for  higher  work  intensities  (M=400  W*m'2, 
-700-800W)  and  lower  ambient  conditions  (Ta=15°C). 

One  factor  that  was  not  addressed  in  this  study  and  awaits  further  study  is 
the  importance  of  the  relationship  of  sweating  rate  and  VO2  max.  Greenleaf  et  al. 
(17)  found  a  linear  relationship  between  these  factors  during  2-h  exercise  bouts 
in  the  heat.  Also,  compared  to  unconditioned  subjects,  physically  fit  subjects 
(V02  max  >  65  mhkg'^min'1)  showed  some  230%  greater  sweating  rates  than 
the  unconditioned  subjects.  Applying  these  results  to  the  present  study  suggests 
that  highly  fit  subjects  would  saturate  the  thermoregulatory  sweating  response 
system  faster  and  Ereq  to  achieve  heat  balance  would  occur  earlier.  However, 
during  uncompensable  heat  stress  conditions,  excessive  sweating  with  little 
evaporation  would  lead  to  early  water  saturation  of  a  given  clothing  system, 
possibly  changing  the  latent  heat  loss  dynamics  through  the  system.  It  is 
uncertain  whether  the  above  prediction  equations  would  be  supported. 

Another  factor  not  addressed  by  this  study  was  the  sweat  prediction 
analysis  using  another  model’s  output,  briefly  explained  in  the  Introduction.  The 
SCENARIO  model  uses  the  thermoregulatory  control  sweating  equation  put  forth 
by  the  work  of  Nadel  and  co-workers  (24).  The  equation  is  embedded  in  the 
model’s  code,  and  sweating  rates  are  solved  based  on  changes  in  core 
temperature  and  skin  temperature  output  using  the  following  format,  where  the 
prediction  of  sweating  rate  (g*min_1)  is: 

msw  =  AD*c))sw,[ot(Tbi  —  Tbi,0)  +  P(TSk  — Tsk!o)]®xp  (TSk_TSk,o)/io 

(Eq.  6) 

In  the  equation  above,  Ad  is  BSA  area  (m2),  a  is  the  control  coefficient 
(4.83  gTnin'^C'1)  modifying  the  change  in  blood  temperature  (i.e. ,  bathing  the 
hypothalamic  thermoregulatory  center,  Tbl,x,  °C)  at  a  given  time  based  on 
exercise  or  other  factors  from  an  arbitrary  hypothalamic  set  point  (Tbl,o  = 
36.96°C),  and  p  is  the  control  coefficient  (0.56  gTnin'^C'1)  driving  the  change  in 
average  skin  temperature  (Tsk,x,  °C)  from  a  thermoneutral  set  point  (Tsk= 
33.0°C). 

The  parameter  c))sw  (units  of  gTn'^min'1)  in  the  SCENARIO  model  is  a  new 
critical  factor  driving  all  variables.  This  factor  was  originally  based  on  the  linear 
relationship  between  sweating  rate  output  (gTn'^min'1)  induced  by  pilocarpine 
iontophoresis  and  V02iriax  (2).  The  relationship  was  evaluated  for  both  men  and 
women  with  maximal  aerobic  capacities  ranging  from  33  to  78  mbkg'^min'1.  In 
this  format,  we  see  that  (j)sw  =  0.16*VO2max  -3.16  (and  the  maximum  aerobic 
power  of  an  individual)  modifies  explicitly  the  overall  sweating  output  determined 
by  the  changes  in  internal  body  and  skin  temperatures.  The  above  equation  has 
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only  been  validated  for  minimally  dressed  individuals  and  may  not  hold  for  wider 
environmental  challenges.  More  studies  are  called  for  investigating  this  crucial 
property  on  the  control  of  thermoregulatory  sweating. 

Indeed,  by  using  the  wealth  of  data  on  2-mile  run  times  (generally  used  in 
U.S.  Army  physical  fitness  testing)  garnered  by  the  U.S.  Army  (19),  it  would  be  a 
straightforward  undertaking  to  predict  V02imax  (ml»kg"1*min"1)  using  the  following: 

for  women:  V02max  =  72.9  -  1 ,77x  (r=  -  0.892)  (Eq.  7) 

and  for  men:  V02max  =  99.7-3.35x  (r=-0.906)  (Eq.  8) 

where  x  is  the  given  2-mile  run  time  completion  (min). 

The  parameter  c|)sw  could  then  be  determined  from  these  data  and  used  as 
input  in  the  above  algorithm  to  derive  sweating  rate  and,  therefore,  the  required 
fluid  intake  at  V02max,  or  a  percentage  of  the  aerobic  capacity  to  accomplish  a 
given  military  task. 

Montain  (22)  found  that  the  sweating  output  using  the  SCENARIO  model 
also  overpredicted  responses  compared  to  measured  sweating  rates,  averaging 
some  0.27,  0.58,  and  0.57  L/h  less  than  predicted  for  WBGT  environments  of 
25.6°,  29.5°,  and  33.3°C,  designated  as  heat  category  I,  III,  and  V  conditions, 
respectively.  Additionally,  in  comparing  the  HSDA  and  SCENARIO  models,  Xu 
(unpublished  report,  2007)  found  model  prediction  levels  were  some  20%  higher 
than  observed  values  using  averaged  data  from  Cheuvront’s  study  (4). 

All  the  above  results  would  suggest  that  output  from  prediction  models 
should  be  employed  with  caution  and  not  hastily  used  outside  the  limited  range  of 
climatic  conditions,  clothing  variables,  and  physiologic  and  heat  exchange 
variables  that  were  not  originally  specified  when  formulating  a  mathematical 
model. 
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CONCLUSIONS 


•  In  conclusion,  the  results  from  the  current  project  verify  that  the  original 
equation  used  to  forecast  sweat  loss  over  wide  heat  stress  conditions 
predicts  too-high  values  for  the  data  set  evaluated,  as  shown  in  Figures  5- 
6. 

•  The  original  sweat  loss  equation  (OSE)  =  27.9*Ereq»(Emax’455)  is  based  on 
a  limited  data  set,  old  static  manikin  analyses,  restricted  to  male 
responses,  and  short  exposure  times  (<  120  min). 

•  Caution  should  be  taken  when  using  any  prediction  models  to  develop 
guidelines  for  fluid  replacement  based  on  sweat  loss  prediction. 

•  The  tactic  used  in  a  fluid  replacement  recommendations  report  by  Montain 
et  al.  (22)  is  a  good  rational  approach,  whereby  a  conservative 
assessment  is  taken  to  prevent  over-hydration  by  not  drinking  more  than 
12  L*day_1  coupled  to  a  given  heat  stress  condition  and  work  intensity  of 
the  individual. 


RECOMMENDATIONS 

It  is  recommended  that  the  findings  within  this  report  be  used  to  substitute 
or  modify  the  current  HSDA  program  with  corrected  algorithm  s  below: 

Sweat  loss  (g*m'2*h'1)  =  147+1 .527»(Ereq)  -  0.87»(Emax)  (Eq.  2) 

This  equation  was  found  to  take  into  account  effects  of  heavy  work  and 
clothing  factors,  body  armor  effects,  longer  exposure  times  (2-8  h),  residual 
errors  inherent  in  the  original  equation  to  predict  sweat  loss,  and  is  not  gender 
specific. 

Alternatively,  a  correction  to  the  original  equation  can  be  used  as  a  simple 
replacement,  determined  as: 

Sweat  loss  =  147*exp  (0.001 2*OSE)  (Eq.  5) 

These  equations  require  testing  over  wider  effective  radiant  loads  in  the 
field  (effect  of  Solar  Load,  Shapiro  et  al.,  (32))  using  a  larger  database.  The 
findings  within  this  report  serve  as  a  first  order  approach  that  may  allow  better 
prediction  of  individual  water  requirements  and/or  sweat  loss  response  useful  in 
field  trials. 
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APPENDIX  A 


The  heat  balance  in  W*m'2  of  BSA  (AD,  m2)  is  expressed  by  (9): 

S  =  M  -  Wk  -  E  -  (R  +  C)  (Eq.A-1) 

where, 

1 S  =  the  time  rate  of  change  of  body  heat  (gain  or  loss); 

M  =  the  rate  of  metabolic  heat  production; 

Wk  =  the  rate  of  accomplished  mechanical  work; 

E  =  the  rate  of  evaporative  heat  loss  via  regulatory  sweating  from  eccrine  sweat 
glands,  diffusion  (Edif),  respiration  (Eres),  and  metabolic  heat  loss  mr; 

C  =  the  rate  of  convective  heat  loss  from  the  total  body  surface  and  respiration; 

R  =  the  rate  of  radiant  heat  loss  (or  gain  from)  the  surrounding  surfaces. 

RADIATION  EXCHANGE 

In  any  thermal  environment,  a  linear  radiation  transfer  coefficient  may  be 
derived  (9,25)  by: 

hr  =  4-a-a-(Ar/AD)-fC|-  (5.67  x  10  ’8)  [(T0  +  Tsurf)  /  2  +  273.15]3, 

(Eq.  A-2) 


where, 

a  is  the  skin  or  clothing  absorptance  for  the  radiation  exchange  to  the  ambient; 

a  is  the  Stefan-Boltzmann  constant  (5.67  x  10  '8,  W-nT2-K'4); 

the  factor  Ar/Ao  represents  the  ratio  of  the  effective  radiating  area  of  the  human 
body  to  the  total  BSA,  as  measured  by  the  Dubois  surface  area  formula  (=0.72 
for  standing  individuals).  The  interior  environmental  temperature  is  composed  of 
an  average  of  the  operative  temperature  +  all  the  surface  temperatures  (TSUrf) 
including  any  clothing  surface  temperature  (Td);  fd  represents  a  factor  which 
increases  the  effective  Arof  the  body  surface  by  some  15%  per  clo  (9,18).  Mean 


^he  heat  balance  equation  is,  by  definition,  expressed  as  rate  of  exchange  overtime,  so  the 
italics  replace  use  of  an  overdot  in  this  report. 


34 


clothing  surface  is  derived  by  Tsk  in  the  relationship  [T0  +  Fcr(Tsk-T0)]-  For 
shorts+T-shirt  Fci=1 ,  and  Td  =  fsk . 

CONVECTIVE  EXCHANGE 

This  is  represented  by  free  convection  and  forced  convection  via 
increased  metabolic  activity,  or  increased  room  air  movement  artificially.  Two 
equations  for  estimating  the  convective  heat  transfer  coefficient  have  been 
formulated  based  on  a  composite  of  free  and  forced  convection  (25): 

hc  =  1 .2  [  (M  -  50)  (  Pb/760)]  0  39,  W-m^-X'1  (Eq.  A-3) 

where,  M  is  the  metabolic  activity  in  W-m'2  and  PB  is  the  barometric  pressure  in 
Torr  (1  kPa/7.5Torr).  Alternatively,  hc for  fan  generated  forced  convection,  in 
which  ambient  air  movement  (V,  m-s'1)  is  the  main  factor  affecting  convective 
heat  exchange,  can  be  expressed  by  either  (hc  in  W-m^-X'1): 


hc  =  8.6  [  V-Pb/760  ] 0  53 

(Eq.  A-4) 

when  persons  are  dressed  in  shorts  and  T-shirts  or  by: 

hc  =  12.7  [V-PB/760]  050 

(Eq.  A-4') 

when  persons  are  clothed  (9,18,25). 
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APPENDIX  B 


Woodcock  and  Breckenridge’s  (1)  procedures  applicable  to  clothing  heat 
transfer  were  used  to  calculate  environmental  heat  exchange.  These  methods 
consider  the  skin,  clothing,  and  environment  as  a  total  system  and  the  constants 
defining  insulation  and  water  vapor  transfer  as  functions  of  effective  air 
movement  (veff).  veff  is  the  sum  of  air  motion  around  a  stationary  object  plus  the 
speed  at  which  the  object  is  moving.  IT  (vefr)  is  the  total  resistance  to  heat  flow  by 
radiation  and  convection  (in  clo  units,  1  clo  is  equivalent  to  0.155  m2K*W'1or 
thermal  conductance  of  6.45  W  m^K'1),  and  im  (veff)  is  the  relative  total  resistance 
to  evaporative  heat  transfer  (zero  to  one,  dimensionless).  In  heat  balance 
calculations,  lm  is  not  used  alone  but  as  a  latent  heat  transfer  coefficient  (im  /  lT), 
evaluated  in  an  articulated,  moving,  sweating  manikin;  this  latter  quantity  is 
considered  as  a  key  dynamic  constant  incorporating  both  heat  and  mass  transfer 
via  “pumping”  through  cuffs,  vents,  and  walking,  as  well  as  relative  permeation 
from  skin  to  ambient,  important  in  total  latent  heat  transfer  efficiency  of  military 
clothing. 

Values  for  It  and  (im  / 1  t)  as  a  function  of  veff  should  always  be  taken  into 
account  in  calculation  of  Ereq  in  the  summed  heat  balance  whenever  there  is  a 
change  in  the  clothing  system  via  exercise  or  use  of  body  armor  (dynamic 
effects).  These  clothing  parameters  can  be  calculated  from  the  following  power 
curves  automatically  estimated  on  the  sweating,  articulated  manikin  used  to 
evaluate  clothing  ensembles  (12): 

lT  =  A*veff  B  and  (im  /  lT)  =  Oveff  D,  (Eqs.  B-1,  B-2) 

where  the  coefficients  A  and  C  are  the  values  for  lj  and  (im/ 1  t)  when  veff  =  1 .0 
nrs'1  and  the  coefficients  B  and  D  are  slopes  of  plots  of  In  (lT)  and  In  (im/  lT)  vs. 

In  (v). 


The  intrinsic  thermal  insulation  value,  W  (v),  is  obtained  by  subtracting  the 
value  of  the  insulation  of  the  air  boundary  layer,  laci,  from  lT: 


lint  =  It  - 


_ 1 _ 

(facl)*  (0.61  +  1  ,87VVeff) 


(Eq.  B-3) 


where  faci  in  Eq.  B-3  is  the  increase  in  surface  area  due  to  clothing  that  is 
estimated  (1 )  using: 


fad  =  1  +  (0.2  •  Ad).  (Eq.  B-4) 

Breckenridge  (1)  defined  the  algebraic  sum  of  the  total  (DRY)  heat  loss  by 
radiant  and  convective  heat  exchange  (R+C),  watts  as: 
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(Eq.  B-5) 


Dry  =  (R  +  C)  = 


6.45*  ADfacI 


lint 


0.61(T sk-T mr)  + Vv  *(Tsk-Ta) 
—  +  facl(0.61  +  1.87Vv) 


INSENSIBLE  HEAT  LOSS 


E  is  determined  by  the  rate  of  sweat  secretion  (msw)  and  the  maximal  rate 
of  evaporative  heat  loss  from  a  fully  wetted  skin  surface  (Emax).  Emax  is  a  function 
of  the  vapor  pressure  gradient  between  the  fully  wetted  skin  surface  and  the  air 
(Ps.sk  -  Pw),  the  evaporative  heat  transfer  coefficient  (he)  and  im,  Woodcock's 
dimensionless  factor  for  permeability  of  water  vapor  through  clothing.  The 
evaporative  heat  transfer  coefficient,  he,  is  directly  related  to  the  convective  heat 
transfer  coefficient,  hc,  by  the  Lewis  Relationship  (LR,  2.2°C/Torr  or  16.5  K/Kpa) 
(9). 


When  evaporation  is  not  restricted  by  clothing  or  the  environment,  then: 


Esk  =  itlsw  •  A  (Eq.  B-6) 

where  rhsw  is  in  g*h_1  and  A  is  the  heat  of  vaporization  for  sweat  at  35°C  (0.68 
W*h*g_1),  Wenger  (37).  The  expression  for  ESk  under  conditions  where 
evaporation  of  sweat  is  restricted  and  there  is  frank  dripping  (Edrip)  or  wasted 
sweat  due  to  skin  wettedness  (co)  >  1 .0  is  (8): 


Esk  =  (0.06+0.94*  co)  •  Ad  EmaX 
when  Emax  ^  m  SW  *  A 


(LR-6.45)  Ad  *(im/lT  )*(Ps,sk-Pw  ), 

(Eq.  B-7) 


where  AD  is  the  DuBois  surface  area  (m2)  (1,9):  Ps,Sk  (in  Torr)  is  the  vapor 
pressure  of  saturated  air  at  skin  temperature.  Ps,Sk  is  related  to  Jsk  by  the 
Antoine  Equation  (9): 


Ps,  sk 


exp  (  (18.6686- 


4030.183  , 
Tsk  +  235 


(Eq.  B-8) 
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RESPIRATORY  HEAT  LOSS 


( Cres+Eres )  is  directly  related  to  ventilation  rate  which,  in  turn,  is  directly 
related  to  aerobic  exercise  intensity  (Mtot)  up  to  maximal  levels.  The  combined 
equation  for  estimating  respiratory  loss  by  convection  and  evaporation  is  taken 

from  (21)  for  high  levels  of  exercise  as  0.019*  Vo2(44-Pw)  in  g/min): 

(Cres  +  Eres)  =  AD  •  Mtot  [0.001 4*  (34  -  Ta)  +  0.0023  (44  -  Pa)]  (Eq.  B-9) 

Eres  is  also  modified  by  a  constant  (F)  for  high  levels  of  exercise  (28),  so  that  if 
Vo2  <  2.6  L/min,  F  =1  and  if  >  2.6,  F  =  1  +0.106(  Vo2-2.6)2.  (Eq.  B-10) 

EVAPORATIVE  HEAT  LOSS  CORRECTION  FOR  METABOLIC  HEAT 
LOSSES  BY  C02  AND  02  ( m  ) 

•  • 

If  the  Respiratory  exchange  ratio  (R)  is  =  1,  m,es  =  Vo2(R*pC02  -p02), 
g/min,  (Eq.  B-11) 

where  pC02  =  density  of  C02=  1 .96  g/L  STPD, 

and  p02  =  density  of  02=  1.43  g/L  STPD  and  mres=  0.53»Vo2;  if  R  >  1, 

in  res  —  V02  (R*0.53)  (Eq.  B-12) 
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